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Abstract 

We study the steady state of an assembly of microtubules in a confined volume, analogous to the 
situation inside a cell where the cell boundary forms a natural barrier to growth. We show that the 
dynamical equations for growing and shrinking microtubules predict the existence of two steady 
states, with either exponentially decaying or exponentially increasing distribution of microtubule 
lengths. We identify the regimes in parameter space corresponding to these steady states. In the 
latter case, the apparent catastrophe frequency near the boundary was found to be significantly 
larger than that in the interior. Both the exponential distribution of lengths and the increase in 
the catastrophe frequency near the cell margin is in excellent agreement with recent experimental 
observations. 
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Microtubules are long, rigid polymers which play an important role in several cellular 
processes. They are highly dynamic structures, which are perpetually in a state of growth 
or shrinkage, and switch stochastically between these states. This behavior is called dynamic 
instability^, |^. The basic monomer unit of a microtubule is a dimer of a and (5 tubulin, 
which is approximately 8 nm in length. The a — /5 dimers are arranged head-to-tail along a 
microtubule in protofilaments (usually 13 per microtubule). 

The highly dynamic nature of microtubules originates from the hydrolysis of /3-tubulin 
bound GTP. Following hydrolysis, the GTP is converted to GDP jsl, and the GDP-bound 
tubulin does not polymerize as well as its GTP counterpart. When the advancing hydrolysis 
front reaches the growing end of a microtubule, it starts de-polymerizing and the protofil- 
aments start peeling off, releasing the GDP-tubulin into solution. Outside a microtubule, 
reverse hydrolysis takes place and polymerization events start all over again. Thus the 
microtubule constantly switches between phases of growth and shrinkage. The stochastic 
transition from growth to shrinkage is called catastrophe and the reverse transition is called 
rescue. The relative rates of catastrophe and rescue, combined with the velocities of growth 
and shrinkage determine the character of a given population of microtubules^, 0, 0]. 

The first theoretical model of microtubule dynamics based on the dynamical instability 
mechanism was developed by Dogterom and LeiblerjTj. In this model, microtubules are 
assumed to nucleate and grow from a flat substrate, and the dynamics is characterized by 
the velocity of growth (vg) and shrinkage (vs), and the frequencies of catastrophe (z/c) and 
rescue (z^a). In the absence of any boundary which restricts the growth, a steady state 
was achieved when unpg < ucPs, characterized by an exponentially decaying distribution of 
lengths. When this condition was not satisfied, no steady state was reached, and the length 
distribution was Gaussian, with mean length increasing linearly with time, and the width 
evolving diffusively. 

Inside cells, the microtubule growth is constrained by the presence of the cell boundary. 
Experimental observations have shown that the parameters of dynamics show a strong de- 



pendence on the proximity to cell boundarylSl. In particular, the catastrophe frequency was 
markedly higher near the periphery, compared to the cell interior. The obvious explanation 
for this difference was that the growing microtubule loses its 'GTP cap' upon hitting the 
cell boundary and is transformed to a shrinking state. In addition, the length distribu- 
tion of microtubules was found to be exponentially increasing, with a possible dip near the 
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boundary. 



In this brief communication, we show that the exponentially increasing length distribu- 
tion of microtubules can be understood from the Dogterom-Leibler equations, and is a new 
steady state which is a direct consequence of the presence of the cell boundary. We compute 
the steady state distribution exactly, and find excellent agreement with experimental ob- 
servations. We also show that the observed increase in the apparent catastrophe frequency 
near the cell margin can be understood quantitatively within this model. 

Let us consider a set of mcrotubules nucleating from a substrate, and growing by the 
addition of tubulin dimers in the direction perpendicular to the plane of the substrate (the z 
axis). For simplicity, we ignore the three-dimensional structure of individual microtubules, 
and treat them as onc-dimcnsional polymers. Nucleation is assumed to take place at empty 
nucleation sites at a rate v. A microtubule in growing state adds T-tbulin at a rate Pg per 
unit time, and a microtubule in shrinking state loses tubulin at a rate per unit time. 
Also, a microtubule in growing state switches to shrinking state at a rate Vc (catastrophe 
frequency) , and a microtubule in shrinking state switches to growing state at a rate vr (rescue 
frequency). Both rescue and catastrophe are assumed to be purely stochastic events. 

Our principal aim in this paper is to study explicitly the steady state of the system in 
the presence of a boundary. We assume that this boundary is located z — I* . We denote 
by p+(/, i) the fraction of sites in the substrate which has a microtubule of length I at time 
t in growing state, and p-{l,t) denote the same fraction in shrinking state. By convention, 
the fraction of vacant sites in the lattice at time t is denoted p_(0, t), and p+{0,t) = at 
all times t. The discrete equations for the dynamics of this assembly, including growth, 
shrinkage, catastrophe and rescue events are given by 



dt 



i/p_(0,t)+p,p_(l,i) ;i = 



(1) 



and 



dt 



up_{0,t) - pgp+{l,t) - ucP+{l,t)i^RP-{l,t) ;/= 1 



(2) 



dp+{l,t) 
dt 



Pg\p+{1 -l,t)- P+{1, t) + iyRP-{l, t) 



iycP+{l,t) ;l<Z<r 



(3) 
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dv (I t) 

^ = Ps[p+{1 + l,t)-p4l,t) + iy,p4l,t) - URp4l,t) ;!</</* (4) 

The presence of the boundary affects the dynamics of the system in the following way: 
When a growing microtubule reaches a length /*, it is instantaneously transformed to the 
shrinking state with length /*. The equations representing this process are give by, 

= p,p+(r - i,t) ~PsP-{i\t) (5) 
p+(r,t) = o 

To find the steady state of the system, let us put all time derivatives to zero. From EqU] 
and EqI21we get the foUwing relations. 



p_(0) = V(l) (6) 
P^{1) = -P-(^) + -^P-i^) (7) 

Pg + J^c 

After combining EqEland EqEl we find that 

Ps + l^R 



and, after using EqE] again, 

p_(0) = (9) 

For / > 1, we find the following relation between p+{l) and P-{1) from EqlS] and Eq0] 
(using only / > 1 in EqH}. 

p.{i + i)-p-{i) = ^[p41)-p4i - 1)] ^ > 2 

Ps 

The general solution of this equation is 

p_(/) = %+(/-l) + C l>2 (10) 

Ps 

where C is an unknown constant. After substituting Eq^Jin Eq„ we obtain the following 
equation for p+(/). 



Ps 

A trial solution to this equation has the form 



(11) 



p+{l) = Aa' + B ;/>l (12) 

We now substitute this solution into Eq llll and after equating terms with the same power 
of I, we obtain the following expressions for the constants a and B. 



1 + ^ 

Ps 



1 + ^ 

Pg 



I^RPs 



-c 



B = 

J^cPs - ^RPg 

The constant C may now be determined as follows: From EqEl for / = 1, we have 



(13) 
(14) 



Ps[p-i2) -p-il) = urp_{1) + ucp+il), 
whereas from Eq^Jwe have another relation: 



(15) 



p_(2) = ^p+(l) + C 

Ps 



(16) 



We now substitute Eq|H]and EqElinto Eqll5l and solve for C, which gives C = 0. From 
Eq ]141 this also implies B = 0. It remains to determine the constant A, which is found using 
normalization: 



(17) 



which may be written as 



p_(o)+P-(i) + (i + ^)ClV+(0 = i 

Ps 

after using EqllOl We now use EqlHl EqlHl and EqfT^ in Eq^l The final result is 



(18) 



A 



Pa 



+ 



Ps + C-^)[^' 



(19) 



The solution in Eq^J can also be written as 



p+(/)=Ae"' ■,a = log 



The complete length distribution may now be written explicitly: 



1 + 






Ps 


1 + 




P3 . 



(20) 



p(/ = l) = A(^ + a) (21) 

Ps 

p{l) = A{1 + ^)e'^' ;!</</* 
aps 

p{l = /*) = Pl^e'^i* (22) 
Ps a 

If — < — , a < and we have an exponentially decaying solution. On the other hand, if 

Ps Pg 

^ > > 1 and a > 0, and we have an exponentially increasing steady state distribution 
of lengths. 

It is interesting to look at the behavior of the solution in the limit /* oo. From Eq lT^ 
we see that, in this limit, a steady state is possible only if a < 1. For, if a > 1, then A ~ a~^* 
for large I* and vanishes as /* ^ oo. For a < 1 and Z* — * oo, ~° -r—, and so 

' a— 1 1— a' 

Pi>^ y Ps ) l-a 

The exponentially decaying steady state length distribution when /* = oo, with ^ < 

QPs Pg 
. The novel feature in the 

finite /* case is the steady state with exponentially increasing distribution of lengths when 

Ps Pg ' 

The exponentially increasing distribution of microtubule lengths has indeed been observed 
in experiments with real cells. Direct observation of microtubules inside cells has been made 
possible recently 0]. These experiments, done on centrosome-containing cytoplasts, observed 
almost persistent growth of microtubules almost up to the cell boundary. However, the 
catastrophe rate showed a dramatic increase within a zone about 3/im near the cell margin 
(0.08s~^, compared to 0.005 in the cell interior). The other parameters describing the 
microtubule dynamics were ~ 0.12s~^,fg ~ 17.8±13.8/im/mm, ~ 28.8±14.1/xm/mm. 
The parameters Pg and ps are related to Vg and Vs a.s Pg = ^ ^^"^ Ps = where 6 is the 
unit of length for our effective one-dimensional polymers. Since a microtubule consists of 
13 protofilaments, and the length of a single tubulin dimer is 8 nm, this length scale is 
6 = 8nm/13 ~ O.Gnm. 



As a first test of our model, we compute tlie increase in catastroplie frequency near the 
cell margin. Let us consider all microtubules with length between li and I*. The total 
number of such microtubules is given hj N = T!i^p{1). Using the expression for p{l) from 
Eq|221 we find that 

a aps 

The number of microtubules in this set undergoing catastrophe per unit time is given by 

N* = iy,N + pgp4l* -1), 

where the first term is the standard catastrophe term, and the second term represents 
the additional catastrophe events arising from the microtubules hitting the boundary. The 
apparent catastrophe frequency is given by z/* = After substituting for p{l) and p+{l* — 
1), we find that 



jy* = jy^+ ^4 ^ (24) 

aps 

where A = After substituting for all the numerical values and for /* — /i ~ 3/im 

as in experiments, we find that v* ~ 0.0964s'^. This is in excellent agreement with the 
experimentally measured value of 0.08 s~^. 

It is also interesting to compare the experimentally measured value of a with the theo- 
retical value. The observed steady state length distribution was found to fit well with an 
exponential function P{1) ~ e^' with 7^^ ~ 5.8;um[8j. We can convert this value to dimen- 
sionless units by multiplying with our unit of length, which gives a^xp = ^7"^ — 1.03 x 10~^. 
The theoretical value is found from Eql2()| using the measured values of all the parame- 
ters, and turns out to be a 1.5 x 10"^. This is also in in very good agreement with the 
experimental value. The discrepancy between the computed and observed values may be 
attributed to the significant experimental error in the measurements of Vg and Vs- 

To conclude, we have studied the steady state of a microtubule assembly in a confined 
geometry, where the growth of individual microtubules is restricted in length. We found 
that, in addition to the exponentially decaying length distribution in an infinite system, 
there is a novel steady state with exponentially increasing distribution of lengths. This 
prediction is is excellent agreement with experimental observations in real cells, and is thus 
a direct verification of the dynamical instability model of microtubule dynamics. 
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